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1 Introduction. 



The derivation of the thermal properties of a black hole is typically carried 
out in the context of the Euclidean path integral approach, as initiated by 
Gibbons and Hawking |l|]. In this language, which is manifestly geometrical, 
the black hole partition function is identified with the Euclidean path inte- 
gral, the integration over four-metrics playing the role of the configurational 
sum, and the hamiltonian, a function of the metric and curvature, taken to 
be that of the gravitational system. 

A black hole of mass M has a temperature Too = h/{8nM) measured at 
large distances from the hole. Following and extending this approach, York 
demonstrated that the canonical ensemble with elements of radius r and 
temperature T(r) for hot gravity with black holes is well-defined 0. That is, 
one treats a collection of spherical cavities of radius r with a temperature T at 
r. These cavities may contain either, no black hole, or one of two physically 
distinct black holes depending on the value of the product rT. In the case 
when the two distinct solutions pertains, only one of them will correspond 
to a thermodynamically stable black hole. This ensemble resolves a number 
of difficulties in assessing the physical significance of the classical black hole 
action in its contribution to the Euclidean path integral. 

However, the reasons for considering the implantation of a black hole in a 
finite cavity, or "box" , go well beyond the resolution of these initial difficul- 
ties and in fact the spatial cut-off provided by the cavity has been recognized 
as being crucial for making sense of black hole thermodynamics in general, 
quite independent of the path integral approach. For example, when one 
comes to consider the back-reaction of the radiation on the spacetime ge- 
ometry of the black hole, the system comprised of black hole plus radiation 
must be put into a finite cavity, lest the radiation in a spatial region that is 
too large collapses onto the hole, thereby producing a larger one . Related 
to this (but much more general) is the fact that the usual thermodynamic 
limit in which one invokes the infinite system volume limit does not exist 
for equilibrium, self-gravitating systems at a finite temperature. This follows 
since the system is unstable to gravitational collapse, or recollapse, if a black 
hole is already present. This, in practice, presents no problem since physi- 
cally, one only requires that the system can in principle become sufficiently 
large so that fiuctuations become negligible. This peculiarity of gravitational 
thermodynamics will play an important role in the present paper. 
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While the Euchdean path integral approach is well-defined and allows one 
to obtain the same value of the entropy as required by the thermodynamic 
analogy, the Bekenstein-Hawking entropy, it does not shed any light on the so- 
called dynamic origin of the entropy, nor does it explain the "microphysics" 
giving rise to the macroscopic thermodynamic functions, such as the internal 
energy, heat capacity, the equation-of-state, etc., that characterize the black 
hole. This state of affairs has spawned numerous efforts to understand the 
dynamical, or statistical mechanical, origins of black hole thermodynamics 
with particular emphasis paid to an explanation of the dynamical origin of 
black hole entropy as for example in Refs. [0 and 

In contrast to on-going efforts devoted to identifying the black hole's fun- 
damental degrees of freedom, we wish to take a model-oriented approach and 
promote a phenomenological analogy between black holes and liquids. The 
analogy will be established at the level of thermodynamics. In the present 
paper, we seek what might be termed an effective atomic picture of black 
hole thermodynamics. By this we mean that it may be possible to reproduce 
(some) black hole thermodynamics in terms of microscopic properties of an 
interacting fluid or gas. The components of this analogue fluid are massive 
point particles interacting mutually via a pairwise additive potential. If such 
a correspondence is possible, we will have effected a mapping between the 
inherently geometrical degrees of freedom part and parcel of the Euclidean 
approach, and the "atomic" variables actually appearing in standard parti- 
tion functions for fluids. The geometric quantities include metric, curvature, 
manifolds and boundaries. The so-called atomic quantities include the par- 
ticles, their momenta and positions, and their interaction potentials. The 
correspondence is established via the black hole's thermodynamic functions, 
as derived from the standard Euclidean path integral approach, using these 
as given input. The task is then to characterize a liquid or (dense) gas whose 
microscopic properties (as encoded for example by the potential, pressure, 
pair-correlation function) can be "tuned", or adjusted suitably, so as to re- 
produce mathematically the same set of black hole state functions. If the 
program so described is successful, then one has in effect, replaced integra- 
tion over metrics by an integration over a multi-particle classical phase space 
while the gravitational action is replaced by a particle hamiltonian, contain- 
ing a (non-relativistic) kinetic energy and potential energy term. 

In the next Section we review and comment on the essential features of 
black hole thermodynamics as derived from the Euclidean path integral in 
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the saddle point approximation. The black hole energy, entropy, equation 
of state and compressibility are displayed and their qualitative features are 
revealed through various limits and graphical representations. The way in 
which we establish a connection between liquids and black holes is taken 
up in Section III. The key in building a "liquid" model is provided by the 
fundamental equations employed in the study of the atomic dynamics of 
simple liquids and the atomic picture of the thermodynamics of the liquid 
state. These fundamental relations equate the macroscopic (thermodynamic) 
to the microscopic (internal structure, potential energy) properties of fluids. 
These relations are derived as rigorous consequences of statistical mechanics 
applied to fluids and (dense) gases. A particular type of fluid is singled out 
the moment we identify the macroscopic variables of the fluid with those of 
the black hole. The points of contact between the black hole and fluid are set 
up via their respective internal energies and compressibilities. The analog 
fluid is identified to the extent that we can write down its pair-potential 
and two-body correlation function. The, necessarily, bounded spatial extent 
of the black hole ensemble is crucial in allowing us to solve for the liquid's 
microscopic parameters exactly. These are calculated in closed form as well 
as presented graphically. 

The ultimate purpose of establishing such a mapping is the double benefit 
to be gained in being able to relate black hole physics to the molecular 
dynamics of fluids. Recent work of a similar spirit includes the possible 
correspondence between black holes and quantized vortices in superfluids [H 
and a connection between fluid surface tension and black hole entropy [1^ 



A summary is given in Section III. Absolute units {G = c = h = ks = 1) 
are used throughout except where restoration of conventional units may be 
helpful. 



2 Black hole thermodynamics in brief. 

In deriving gravitational thermodynamics from an Euclidean path integral 

Z{(3) = J d4g,4>]e-'^^''^^ = (1) 

one expects the dominant contribution to the canonical partition function to 
come from those metrics g and matter fields that are near a background 
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metric g^^^ and background field 0*-°'', respectively. These background fields 
are obtained from solutions of the classical field equations. The classical 
contribution to Z is obtained by evaluating Z at the point (g^^\(f)^^^) in 
which case one obtains the familiar relation 

/?F = /[^W,0(°)] = J[^W], (2) 

where we have taken 0*^°^ = in the last equality. This provides the free 
energy F of the gravitational system in the saddle-point approximation. The 
action / is the first-order Euclidean Einstein action including a subtraction 
term necessary to avoid runaway solutions. The action appropriate to a black 
hole in a spherical cavity of radius r is given by 

I = h — I subtract, (3) 

where 

I, = --L r f^*d'x^R+^i d'x VlTr{IC) (4) 

The action receives in general both volume and boundary contributions. The 
Euclidean four-space metric is 

9,u = diag ((1 - (1 - ^)-^r^r^in2 0) . (5) 

For this metric, the volume contribution vanishes identically. The boundary 
at r = const, is the product of 5*^ x S"^ of the periodically identified Euclidean 
time with the two-sphere of area A = Anr"^. The period of Euclidean time, 
identified with the S'^-coordinate, is /?* = SvrM. The trace of the boundary 
extrinsic curvature is denoted by Tr(/C) and 7 is the induced 3-metric on 
the boundary. Finally, Subtract is h evaluated on a flat spacetime having the 
same boundary x S^. 

It is important to remember that for the canonical ensemble the mass 
parameter M appearing in these formulae is not simply a constant but is 
instead a specific function of the cavity radius r > and the cavity wall 
temperature T(r) > 0, ^. This can be verified by inverting the expression 
for the blue-shifted temperature in equation (11) below and solving for M = 
M{r,T). The relation so obtained is a cubic equation in M. When rT < 
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\/27/87r = {rT)min ~ 0.207, there are no real solutions of this equation. On 
the other hand, when rT > {rT)min) there exist two real and non-negative 
branches given by 



M2{r,T) 

M,(r,T) 

cos(q;) 




l + 2cos(-) 



= 1 
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(6) 

(7) 

(8) 
(9) 



This shows that the Schwarzschild mass is in fact double- valued in the canoni- 
cal ensemble. One has that M2 > Mi, with equality holding at rT = {rT)min- 
The heavier mass branch, M2, is the thermodynamically stable solution be- 
cause it leads to the lowest free energy, Eq. (2), and is the one we shall be 
considering in the remainder of this work. 

Calculating the action I from (3) and (4) yields {Igubtract — —f^f) 

I = 127rM2 - 87rMr + /?r. 



where 



13 = T- 



SttM 1 



2My/2 



(10) 



(11) 



is the inverse local temperature and is the proper length of the component 
of the boundary. Employing / and the saddle-point approximation (3F = I,it 
is a straightforward exercise to calculate the thermodynamic state functions 
associated with a black hole in the canonical ensemble. In so doing, it is 
useful to note the following two identities 
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(13) 



which may be deduced from the expression for the inverse local temperature 
(11). The black hole's internal, or thermal, energy is 



(14) 
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The entropy S is 



(15) 



while the surface pressure a is 



a 




(i-f) 



^1 _ 2M^ 



1/2 



(16) 



Another quantity of special interest in what is to follow, is the black hole 
isothermal compressibility, kt{A), which again can be calculated using the 
standard prescription. 
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Although at face value these functions appear to have a complicated depen- 
dence on r and T, they are actually quite simple, owing to their dependence 
on the slowly varying ratio M/r. To gain some insight into the behavior 
of these functions, it is useful to examine some limiting cases, namely when 
(i) rT — > oo and when (ii) rT — > {rT)min- The hmit (i) is understood to 
mean either that r — cxd or T — > cxd, or both, simultaneously. The second 
limit {ii) actually defines a hyperbola in the r — T plane along which the two 
independent limits r — >■ (T — oo) or r ^ oo (T — ^ 0) can be taken. The 
mass function takes the form M(r, T) f» |(1 — ^^2^2-7^2 ), and M(r, T) = |, 
respectively. Physically, these hmits indicate that for all allowed values of 
r and T, the cavity wall radius always lies between the black hole's event 
horizon (r = 2M) and the unstable circular photon orbit (r = 3M). The 
behavior of the black hole internal energy with respect to these limits is 



E 



1 



(47rr) ' 



and 



(ii) 



E = r. 

V3 



(18) 



(19) 
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E is essentially a positive linear function of r, depending only very weakly 
on the temperature for large values of rT. For rT = {rT)min, E is strictly 
linear in r, or inversely proportional to T, since of course, in this latter limit, 
r ~ 1/T. Note that on this hyperbola, — > for r — > (or T oo). 
The equation for the surface pressure is also an equation of state for the 
black hole pressure since it is expressed as a function of the cavity radius r, 
which gives a measure of the system size, and the boundary temperature T: 
(7 = (j(r, T). Using the limiting forms of the mass function, one can show 
that the asymptotic limit of the surface pressure is 

for rT — > oo, so that the pressure increases with the temperature, depend- 
ing only very weakly on the system size. When evaluated along the limit 
hyperbola, one obtains 



in other words, in this regime, the pressure increases as the cavity radius 
(or area) decreases. Because of the reciprocal relation between r and T 
along the hyperbola, this is equivalent to increasing pressure with increasing 
temperature. Such qualitative behavior is familiar from the ideal gas. Finally, 
the limiting forms of the isothermal compressibility are 

(i) KTiA) -167rr < 0, (22) 

and 

(ii) kt{A) - 0. (23) 

These two latter limits deserve some special comment. First, note that the 
black hole isothermal compressibility is generally negative. This is an unfa- 
miliar property in regards to conventional thermodynamic systems. Indeed, 
standard textbook arguments prove that kt > 0, irrespective of the nature 
of the substance comprising the system. However, a key step in those proofs 
assumes that quantities such as the temperature and pressure are intensive, 
that is, independent of the size of the system and constant throughout its 
interior. Such is most emphatically not the case for gravitating systems such 
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as black holes, where in fact, the temperature and pressure are not intensive 
quantities but are instead scale dependent. An equilibrium self-gravitating 
object does not have a spatially constant temperature. This is a consequence 
of the principle of equivalence which implies that temperature is red- or blue- 
shifted in the same manner as the frecuency of photons in a gravitational field. 
Secondly, for values rT = {rT)jnin, the compressibility vanishes identically. 
This qualitative behavior is familiar from the classical picture of a solid at 
T = (no density fluctuations =^ zero compressibility). 

In the Figures 1-3, we have plotted E, a and in the r — T plane subject 
to the condition rT > {rT)min- As indicated in Fig. 1, the black hole energy 
is a positive increasing function in r and is fairly insensitive to changes in the 
temperature for values of T > 0.4. The flat region at zero level corresponds 
to the locus of excluded points satisfying rT < {rT)min and is therefore not to 
be considered as part of the graph as such. The cavity wall surface pressure, 
shown in Fig. 2, is a positive increasing function of T and varies slowly in 
r for r > 0.5. The flat null region represents the same locus of points as 
in the prior graph. Finally, the black hole isothermal compressibility is a 
negative definite function for all rT > {rT)min, decreasing for increasing r 
and relatively constant with respect to changes in the temperature. 

Other functions that may be calculated include the specific heats at con- 
stant area and at constant pressure, respectively, as well as the adiabatic 
compressibility, but these are of no direct interest for the present considera- 
tion. 

Finally, to complete this brief overview of black hole thermodynamics, 
we need to identify the effective spatial dimension of the system. It will 
not have escaped the reader's attention that the above functions have been 
defined and calculated in terms of the cavity wall area A, rather than in terms 
of the cavity volume. Spatial volume is not well defined in the presence of a 
black hole, whereas the area is, and it is the latter which provides the correct 
means for measuring the size of the system |Q. That the wall area is the 
proper extensive variable to use is confirmed by considering the black hole's 
thermodynamic identity. The explicit calculation of the following partial 
derivatives 




(24) 
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and 




I _ M 



^1 _ 2M 



xl/2 



-a. 



proves that 



dE^TdS-adA 



(25) 



(26) 



is an exact differential. In other words, the energy when expressed in terms of 
its proper independent variables, E = E{S, A), is integrable. We remark that 
all the above functions may be considered as functions either of the cavity 
radius or the wall area, since obviously, r = JA/{A7r) and dA = STirdr. 



3 A Liquid Model for Black Holes. 

The starting point for attempting to model black hole thermodynamics in 
terms of liquids is the statistical mechanical treatment of fluids. There, it is 
known how to relate the various macroscopic thermodynamic properties of 
liquids (energy, pressure, temperature) to the internal, microscopic features 
such as the interaction/intermolecular potential and the pair-correlation func- 
tion. This latter function provides a measure of the local structure of the fluid 
or gas. The typical approach to the study of the liquid state starts with (a 
perhaps imperfect) knowledge of the interatomic force law and the measured 
short range order (obtained experimentally via X-ray or neutron scattering 
experiments). One then attempts to infer the macroscopic or thermody- 
namic behavior of the liquid on the basis of this microscopic information. By 
marked contrast, here we shall turn the reasoning around and solve for the 
local "structure" and the "interatomic" potential of an analog (and possi- 
ble fictitious!) fluid from explicit knowledge of black hole thermodynamics. 
Two very important ingredients which will allow this inverted procedure to 
be carried out in closed form are (i) the fact that the thermal ensemble of 
black holes is spatially bounded and (ii) the fact that we know the spatial 
dimension of this ensemble to be c? = 2. 

The model fluid we shall deduce will be described in classical terms. Let 
us use our knowledge of the spatial dimensionality of the black hole ensemble 
at the outset. To make a thermodynamic correspondence between black hole 
and fluid means we seek a two-dimensional fluid whose partition function 
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over the 2A^-dimensional phase space is given by (restoring the dependence 
on h) 

and where the total, nonrelativistic energy of a system of interacting point 
particles of mass m in d = 2 is 

27V 2 

mP^},{r^})=T.^ + U{{r.})■ (28) 

Here, U{{ri}) is the potential energy of the particle system which we assume 
to be pairwise additive: 

U=lj:<Pi\r^-r,\). (29) 

This is always a reasonable assumption provided the fluid constituents have 
no internal structure that couples to the potential. From the theory of liquids, 
it is well known that an equation of state for the isothermal compressibility 
kt and the internal energy E can be calculated in terms of the pair-potential 
and the pair-correlation function g [|l^, to whit (restoring the depen- 
dence on ks), 

(a) : p ksT^T = P f rf'f W) - 1] + 1, (30) 

J system 



and 



(6) : E = NksT + ^-Np f d^i <P{f)g{f), (31) 

/ J system 



where p = ^ is the two-dimensional particle density, T is the fluid tempera- 
ture, and the radial distribution function g is defined via 

P^?(0 = ^(E^(r + r^-r,)), (32) 

the angular brackets denote the average computed using the grand canonical 
ensemble. 

Before we go on to use these relations, we remark that the equation of 
state (a) is exact while the expression for the internal energy (6) makes use of 
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the pairwise summability of the total potential energy. They are valid for any 
single-component, monatomic system in thermodynamic equilibrium (gas, 
liquid or solid) whose energy is expressible in the form (28) with a pairwise 
additive potential (29). Although these expressions are derived primarily for 
their application to the liquid state, they can also be applied to the study of 
solids. The only modification would be that g and (j) depend on the full vector 
coordinate r (magnitude and direction). For liquids, however, the results are 
isotropic so it is enough to write g{r) and (j){r). For the present consideration, 
modelling a liquid which is capable of reproducing certain aspects of black 
hole thermodynamics is carried out once we identify the kt and E in (a) and 
(b) with those of the black hole. 

The idea of representing a black hole at finite temperature as a thermal 
fluid is novel and therefore deserves careful explanation. A black hole is 
but one example of a thermodynamic system having, among other things, 
a well-defined temperature, energy and compressibility. On the other hand, 
any equilibrium many-body system with hamiltonian given in Eq.(28) and 
Eq.(29) has a well-defined energy and compressibility which can be calculated 
in terms of an associated g and 0. When we formally identify the kt and 
E appearing there with those belonging to the black hole, we are simply 
demanding that these particular thermodynamic functions be reproducible 
in terms of the internal variables of a certain classical many-body system. 
This is not to say that these variables actually represent the true degrees 
of freedom of a black hole. The identification must be carried out in a 
consistent way. First, the temperature T appearing in (30) and (31) is the 
uniform cavity wall temperature. Since (a) and (6) are to decribe a liquid, 
the temperature of that liquid must be identified with this temperature: 
Tiiquid = T. Note that the temperature of the liquid is intensive. That is, the 
temperature of the cavity wall of the black hole ensemble is identified with 
the temperature of the bulk fluid. Next, the density p of the fluid is simply 
the number of "atoms" per unit area of fluid. For the black hole, both E 
and Kt depend explicitly on the cavity radius r, refiecting the fact that the 
black hole ensemble is spatially finite. This means that the integrations in 
(a) and (6) are to be carried out over a fluid of bounded spatial extent. The 
integrations over the system are bounded. Since the integrands are functions 
only of r and T, we can therefore write 
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/ r dr d9. 



(33) 



J system 



It is natural to take the length scale of the liquid coincident with that 
of the cavity containing the black hole; any other choice would introduce a 
second, and arbitrary, length scale into the problem. Taking E and kt from 
(14) and (17) as input, the relations (a) and (6) yield two equations in the 
two unknowns g and 0. We can easily solve for these microscopic functions 
in terms of the macroscopic functions and their first derivatives. To do so, 
we make use of (33) and differentiate both sides of the relations in (a) and 
(6) with respect to r. The results of this operation are that 



By explicit construction, these give the pair correlation function and the 
inter-particle potential of the model fluid whose energy and isothermal com- 
pressibility are identical with those of the black hole. Moreover, these two 
functions depend on the two independent variables r and T. Since we can 
vary them independently, we have actually obtained (p and g as functions of 
their arguments for all T > and r > 0, subject only to the restraint that the 
product rT always be greater than or equal to {rT)min- By way of a trivial but 
illustrative example, consider the ideal gas in two-dimensions whose equation 
of state is pA = NksT. Then E = NksT and kt = = l/{pkBT). 

Inserting these into the above relations immediately yields g{r) — 1 and 
(f){r) = 0, which is also a solution of the pair of equations (30) and (31). As 
is to be expected, the ideal gas has no structure (it is uniform: homogeneous 
and isotropic) and lacks interatomic interactions (by deflnition). Therefore, 
any deviation in either g and or (p with respect to these limits may be con- 
sidered as deviations from an ideal gas. 

It is of interest to consider the limiting forms of the pair correlation 
function and potential energy for the black hole; these may be deduced easily 
from the associated limits calculated above for E and kt, Eqs. (18,19) and 




(34) 



and 




(35) 
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Eqs.(22,23). When rT » {rT)min, the pair correlation function goes as 

W g(r,T)~ 1-^-1 (36) 

In particular, for fixed temperature, g{r, T)— >1 — 5, asr— >oo where S — 
4/N is small for large. In normal simple liquids, g{r) has the asymptotic 
limit g{r) ^ 1 (compare to the ideal gas limit) and deviations from this value 
represent molecular correlations (or anti-correlations). When evaluated along 
the boundary hyperbola rT — {rT)^^^ we get, 

(ii) ^(r,r) = l-l, (37) 

a constant independent of r and T. The corresponding limits for the two- 
body potential energy may be worked out and yield 

(i) 0(r,T)^ 
For fixed temperature, 

(ii) (t>{T,T) 

The black hole pair correlation function is calculated and presented in Fig. 

4. For fixed T and small r, this function is negative, then increases, becoming 
positive and approaches unity from below as r — oo. This behavior is also 
revealed in the one dimensional plot of g{r, T) for the value T — 0.5 in Fig. 

5. What can we make of this behavior in g and what physical interpretation 
can it admit? For this, let us turn to the meaning of g{r). Imagine we 
select a particular particle of the fiuid, whose average density is p, and fix 
our origin at that point. Then, the number of particles dN contained within 
the (two-dimensional) spherical shell of thickness dr centered at r = is 

dN ^ p g{r)2nr dr. (40) 

Here we see that g gives a measure of the deviation from perfect homogene- 
ity {g ^ 1). Evidently, g < <S=^ dN < in that shell. On the other 
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hand, a negative value for dN is the signature for the phenomenon of charge- 
screening, i.e., it indicates the presence of holes in the neighborhood of our 
reference particle at r = 0. Thus, it would appear that the analogue fluid 
which could model some aspects of black hole thermodynamics should have 
something to do with a charged fluid or plasma. These latter systems are de- 
fined as a collection of identical point charges, of equal mass m and charge e, 
embedded in a uniform background of charge (the dielectric) obeying classical 
statistical mechanics. If one adds a given charge to the plasma, the plasma 
density is locally depleted so as to neutralize the impurity charge. This is 
the well-known phenomenon of charge screening. The depletion shows up 
as an underdensity of particles (or an overdensity of holes), and is reflected 
in a 5^ < near the origin, that is, where the impurity charge is located. 
Calculations of the pair-correlation function for degenerate electron plasmas 
at metallic densities yield functions exhibiting the same general quantitative 



features as those in Figures 4 and 5 [|T^]. In addition, screening is known to 
be a characteristic property of interactions like the electromagnetic interac- 
tion, where there exist two species of charges; they have the property that 
renormalization effects induce corrections that make the effective charge 
decrease with distance. This, in itself, is not surprising because as is well 
known there exists an analog of a black-hole with an electromagnetic mem- 
brane; an analogy that in the literature is called the "membrane paradigm" 



IJ] . What we find here is a different manifestation of this analogy, this time 
through the dynamical and statistical properties of the liquid. 

The weakly temperature dependent potential is scaled by A^^ and plotted 
in Figure 6. Again, recall the physical part of the graph consists of those 
points r and T satisfying rT > {RT)min- The potential is seen to be a 
positive increasing function of r. The limit calculated above in Eq.(38) shows 
the growth is essentially linear. Apart from the "glitch" near T ^ 0.2 the 
potential is practically independent of the temperature. 



4 Conclusions. 

It is worth emphasizing that the analogue fluid selected to account for the 
black hole compressibility and internal energy was "engineered" at the fluid's 
atomic level. As there is no corresponding "atomic" level for the black hole, 
the bridge between the thermodynamics of the black hole and liquid is es- 
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tablished via thermodynamic state functions. Surprisingly, only two state 
functions are needed in order to specify completely the "atomic potential" 
and local structure of the analogue fluid. However, as we have seen, we can 
only be sure that the fluid will reproduce the correct compressibility and 
internal energy. That is, only partial aspects of black hole thermodynamics 
will be reproducible, since, evidently, there exist other state functions that 
characterize a black hole, namely, its entropy, pressure, speciflc heats, etc. 
which must be calculable by a more complete "microscopic" description of 
the black hole (see below). The analogy with the liquid leads to a screening 
effect that can be understood in terms of a connection with the membrane 
paradigm. 

The partial rendering of black hole thermodynamics in terms of atomic 
fluid elements achieved here points to the possibility of directly effecting 
a mapping between the black hole variables (mass M and cavity radius r, 
or cavity wall temperature T{r)) and the internal variables of an analogue 
model which might serve to reproduce all of black hole thermodynamics. 
Evidently, this would amount to a formal correspondence at the level of the 
degrees-of-freedom and so bypass the need to call into play the macroscopic 
state functions. A concrete example of such a mapping between two entirely 
distinct systems is that established recently between a Newtonian cosmology 
of pointlike galaxies and spins in a three-dimensional Ising model [0. The 



degree-of-freedom mapping problem is well worth pursuing as intriguing deep 
connections between gravitation, thermodynamics and information theory 



have been hinted at recently [16|. Another hint is supplied by Wheeler's 
depiction of the Bekenstein bit number as a "covering" of the event horizon 
by a binary string code representing the information contained in the black 
hole |18|. It may well be possible to go beyond these provocative hints and 
actually establish a rigorous connection between black holes, computation, 
information theory and complexity. We hope to report on these developments 
in a separate paper. 
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Figure Captions. 



Figure 1. The black hole energy, as a function of the cavity radius r 
and cavity wall temperature T. It is seen that E grows with cavity radius 
and with temperature. The flat portion of the graph (in this and in the 
following figures) is not part of the physical quantity shown in the graph; it 
corresponds to the region rT < {rT)min in which the black hole mass is no 
longer real and positive. 

Figure 2. The black hole surface pressure, cr, as a function of the cavity 
radius r and cavity wall temperature T. It is seen that a grows sharply at 
first with cavity radius to then stabilize its growth, while its growth with 
temperature is unimpeded. The fiat portion of the graph (as in all figures in 
this paper figures) is not part of the physical quantity shown in the graph; it 
corresponds to the region rT < {rT)^^^^ in which the black hole mass is no 
longer real and positive. 

Figure 3. The black hole isothermal compressibihty, kt, as a function of 
the cavity radius r and cavity wall temperature T. It is seen that kt decreases 
with cavity radius and with temperature. The flat portion of the graph (as 
in all figures in this paper figures) is not part of the physical quantity shown 
in the graph; it corresponds to the region rT < {rT)min in which the black 
hole mass is no longer real and positive. 

Figure 4. The black hole pair correlation function, g, inferred in the liquid 
model of the black hole as a function of the cavity radius r and cavity wall 
temperature T. It is seen that for r > .2, gf grows with cavity radius but 
decreases with temperature. The fiat portion of the graph (in this and in the 
following figures) is not part of the physical quantity shown in the graph; it 
corresponds to the region rT < {rT)mi,i in which the black hole mass is no 
longer real and positive. The spikes are numerical artifacts of the plotting 
routine used to plot the figure and which could not render the deep trough 
from which the physical region emerges. 

Figure 5. Two dimensional representation of the black hole pair correla- 
tion function, gi, as a function of the cavity radius r but for a fixed value of 
the cavity wall temperature T — 0.5. This is plotted here to illustrate the 
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existence of the trough which was not clearly seen in the three dimensional 
plot of Figure 4. Notice that for large r the pair correlation function tends 
to the constant value +1. 

Figure 6. The black hole "liquid" two-particle potential, 0, as a function 
of the cavity radius r and cavity wall temperature T. It is seen that apart 
from the sudden rise near T = 0.2, the potential is a rather boring function of 
cavity radius and temperature. The flat portion of the graph (in this and in 
the previous figures) is not part of the physical quantity shown in the graph; 
it corresponds to the region rT < {rT)jnj_n in which the black hole mass is no 
longer real and positive. 
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